Wave functions and decay constants of B and D mesons 
in the relativistic potential model 
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With the decay constants of D and Ds mesons measured in experiment recently, we revisit the 
study of the bound states of quark and antiquark in B and D mesons in the relativistic potential 
model. The relativistic bound state wave equation is solved numerically. The masses, decay con- 
stants and wave functions of B and D mesons are obtained. Both the masses and decay constants 
obtained here can be consistent with the experimental data. The wave functions can be used in the 
study of B and D meson decays. 
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I Introduction 

The wave function of the bound state of quark and an- 
tiquark is determined by the strong interaction between 
quark and antiquarks. The study of the wave functions 
of heavy-flavored mesons like B and D are important not 
only for studying the property of strong interaction be- 
tween heavy and light quarks, but also for investigating 
the mechanism of heavy meson decays. The wave func- 
tion determines the momentum distributions of the quark 
and antiquark in mesons, which is an important quan- 
tity for calculating the amplitude of heavy meson decays 
[11,13] . The light-cone momentum distribution amplitudes 
for B meson appear in the amplitude of B decays, which 
are defined through the hadron-to-vacuum matrix ele- 
ment of non-local operators of quark and antiquark sepa- 
rated along the light-cone {0\q'^{z)b°'{0)\B{p)). The light- 
cone distribution amplitudes of B meson have been ex- 
tensively studied in the recent several years. Some prop- 
erties of the light-cone distribution amplitudes have been 
obtained. Based on these achievements, several models 
satisfying these constraints have been proposed in the 
literature 0-01 • Methods to obtain the light-cone distri- 
bution amplitude exactly from the first principle of QCD 
are still under investigation. 

Alternatively, directly studying the wave functions of 
the heavy mesons by solving the bound state wave equa- 
tion is an effective way to obtain the knowledge about the 
bound state of quark and antiquark [lol - [l5| . For a heavy- 
light system, in the heavy quark limit, the heavy quark 
can be viewed as a static color source in the rest frame of 
the heavy meson. The light antiquark is bound around 
the heavy quark by an effective potential. The heavy 
quark spin decouples from the interactions as uiq — > oo 



{niQ is the mass of the heavy quark) [16|, [17| . The inter- 
actions relevant to quark spin can be treated as pertur- 
bative correction. 

Inspired by asymptotic freedom at short distance in 
QCD and quark confinement at long distance, the effec- 
tive potential between quark and antiquark in meson can 
be taken as a combination of a Coulomb term and a linear 



confining term. Such a potential will be consistent with 
perturbative QCD at short distance, it can also generate 



quark confinement at long distance [ll|, |18| . 



The parameters in the effective potential can be con- 
strained by comparing the eigenvalues of the bound state 
wave equation with the masses of the relevant bound 
states measured in experiment. Recently the decay con- 
stants of D and Dg mesons have been measured by CLEO 
;i9-21]. Belle ^ and BABAR Q CoUaborations. The 
measured values of D and Ds mesons' decay constants 
Jd and Jd^ can give further information about the inter- 
actions within the heavy-light quark-antiquark system. 
In this paper, with the recently measured decay con- 
stants fo and available, we revisit the study of the 
bound states of iJ, Bg, D and mesons in the relativis- 
tic potential model. We solve the relativistic version of 
Schrodinger equation for the bound state wave function 
of heavy- light quark-antiquark meson system. The decay 
constants fo and can be used as a further constriant 
on the parameters in the potential model. The obtained 
masses and decay constants of B and D mesons can be 
well consistent with the values measured in experiments. 
Then the wave functions obtained here can be more re- 
liable then ever. It can be useful for studying B and D 
decays, where the momentum-distribution of the quarks 
is needed. 

Although the bound states of heavy mesons have been 
studied with the relativistic potential model in the lit- 
erature several years before, these works need to be im- 
proved with the recent experimental data of the decay 
constants of D and Dg mesons available. The decay con- 
stants and bound state masses are calculated in Refs. 
[Tsj l , where Richardson potential j23| was taken, here the 
potential we considered is different from theirs. In ad- 
dition, with the experimental values of the decay con- 
stants fu and available recently, the parameters in 
the potential can be constrained more stringently. There- 
fore our prediction on the decay constants for B and D 
mesons are quite different from previous predictions in 
the relativistic potential model. 

The paper is organized as follows. In section II, we 
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solve the relativistic wave equation for the heavy-hght 
quark-antiquark system. Section III gives the decay con- 
stant in terms of the wave function. In section IV the 
QCD-inspired potential is presented. Section V is de- 
voted to the numerical result and discussion. Section VI 
is a brief summary. 



II The relativistic wave equation for heavy-light 
system and the solvement 

The B and D mesons are assumed to be approximately 
described in terms of heavy-light valence-quark configu- 
rations in the rest-frame of the mesons. The effective 
potential is one-gluon exchange dominant at short dis- 
tances and a linear confinement at long distances 
The equation describing the bound state wave functions 
is a Schrodinger-type wave equation with relativistic dy- 
namics 1 1 - 151 



condition 



-n^S/l+ml + V{r) 



(1) 



where f — xi ~ X2 is the displacement of the light an- 
tiquark from the heavy quark, and xi and X2 are the 
coordinates of the light and heavy quarks respectively. 
The operators Vi and V2 are the gradient operators rel- 
evant to the coordinates of xi and X2- mi is the mass 
of the light antiquark, and m2 the mass of the heavy 
quark. V{r) is the effective potential of strong interac- 
tion between heavy and light quark-antiquark. In the 
rest frame of the bound state system, the eigenvalues in 
the wave equation will be the masses of the series bound 
states. 

The wave function can be expressed in terms of spec- 
trum integration 



■0(r) 



^ik-{r—r ')/h 



W)- (2) 



Substitute Eq. ^ into Eq. ([T]) , the wave equation becomes 

xe^-^^-^'^'^iPir') ^{E- y(r))V(r). (3) 

The exponential e^^'^^^ can be decomposed in spherical 
harmonics 



, kr 



(4) 



Ln 



where is the spherical Bessel function, Y;„(r) is 

the spherical harmonics, which satisfies the normalization 



(5) 



Using the spherical harmonics decomposition of the ex- 
ponential in Eq.Q, and factorize the wave function into 
the product of two parts: radial and angular wave func- 
tions 



(6) 



then the wave equation of Eq.Q can be transferred to 
be 



dr'r'^i^ 



+ ^k^ + ml )M^)ji{^)^i(r') = E^r)- (7) 

For convenience later, let us define a new reduced radial 
wave function ui{r) by 



uiir) 



(8) 



With this definition, and for the case I — which we are 
interested in this work, Eq.Q becomes 



V(r)uo{r) + — dk dr'Uk^+mj 
■^n- Jo Jo ^ 

+ ^ k'^+ml ) sm(y ) sin( — )uo(?-') = EuQ{r), (9) 

where we have used the explicit expression of the spher- 
ical Bessel function for I — 



io{x) 



sm X 



(10) 



Eq.® is for taking c = 1, if recover the speed of light 
appearing in the formulas, Eq.® should be 



V{r)u,{r) + — 



dk 



dr' {\J fc2 + m\ 



+ ^J k'^ + ml) sin( — ) sm{ — )uo{r') = EuQ{r). (11) 

In principle the integration over momentum k in the 
above equation can be performed because the wave func- 
tion uo{r') does not depend on the momentum. However 
the integration over k will give a singular term for r' ^ r 
in the above equation IJ]. In this work, we will take a 
new step to continue to solve this equation, this method 
can circumvent the appearance of the singular integral 
equation. 

For a bound state of two particles, when the separation 
between them is large enough, the wave function will 
effectively vanish. We assume such a large enough typical 
value for the separation between the heavy quark and 
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the light antiquark is L, then the quark-antiquark in the 
bound state can be approximately treated as if they are 
restricted in a limited space < r < i. In the limited 
space, the Fourier expansion of the reduced wave function 
uo[r) is 



Uo{r) = ^ c„ sin {^r^ , 



(12) 



n=l 



where the expansion coefficients c„ are 



I sin (^"^''^ UQ{r)dr. (13) 



In the limited space, the momentum k should be dis- 
cretized, the integration over k should be replaced by a 
summation, the following substitution should be made in 
the wave equation (fTTj) 



k mr 
he L ' 



dk TT 
he L 



(14) 



With the above replacement, and the integration over 
the distance r' being limited within < r' < L, Eq.((TT]) 
becomes 



V{r)uo{r) 
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dr'i 



' nirhe 
L 




nnhe 



Euo{r). 



sm( — r) sm( — r )UQ(r ) 



(15) 



The above equation can go back to Ea. ([TT|) as L — ?> oo. 
Numerically if the value of L is taken to be large enough, 
the solution of this equation only slightly depends on the 
value of L. For the parameters we take in section V, 
we find that the solution of the wave equation will be 
stationary when L > 5 fm. 

Truncate the series of the Fourier expansion of the wave 
function uo(r) as 



uoir) 



JV 

E 

ri=l 



c„ sm 



(16) 



where TV is a large integer. Substitute this truncated 
expansion into the wave equation (jlSp and simplify it, 
one can finally get the equation about c„ 




m—l 

■- Eo 



(17) 



The above equation is just the eigenstate equation in the 
matrix form. It is not difficult to solve it numerically. 
The eigenvalues are the masses of the series of bound 
states of the heavy-light quark-antiquark system. Once 
the eigen equation is solved, the eigenvectors composed 
of c„ can be substituted into Ea.([TC| to get the reduced 
wave function Uf){r). 

To get the wave function in momentum space, one can 
use the Fourier transform of the wave function '4){r) 



*(fc) 



1 



(27r?ic)3/2 



d^re-'^-^/'"''ip{r}. (18) 



Separate the variable-dependence of the momentum- 
space wave function as 



^{k) = ^i{k)Yi^{e,^). 



(19) 



As in Eq.®, we define the reduced wave function in mo- 
mentum space 



(20) 



Then using Eqs.dJ), ©, ® , dH]) and one can derive 
from Eq.([T8| 



/ O f°° Ur kr 
Vi{k)^{-i)\— dr-M-)ui{r). (21) 
V nne Jq ne Tie 



For the case / = 0, we get 



drsin(|^)wo(r), (22) 
ne 



which gives the momentum distribution of the quark and 
antiquark in the rest frame of the heavy meson. 



Ill The pseudoscalar bound state of heavy-light 
system and the decay constant 

The pseudoscalar meson composed of a heavy quark 
and a light antiquark Qq [Q can be 6 or c quark, q stands 
for u, d or s quark) can be written in the meson rest frame 
as follows 



1 



1 



|P(p = 0)) = -^^y d3fcvI/o(fc)-^[6^+(A:-,tX(-fc,;) 



-6^+(A:,i)d^+(-A:,t)]|0), 



(23) 



where i is the color index. The factor 1/V3 is the nor- 
malization factor for color indices, and l/\/2 the normal- 
ization factor for spin indices. 

The normalization of the meson state is 

(P(pi)|P(p2)) = {2^f2E6\pi -P2), (24) 

where E is the energy of the meson. 
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Substituting Eq. ((23)) into Ea. ([24)) . we can finally get 
the normalization condition of the wave function in mo- 
mentum space 



rf'fc|*o(fc)|' 



{2t:Y2E. 



(25) 



The decay constant of a pseudoscalar is defined by the 
hadron-to- vacuum matrix element of the axial current 



]qit^ibQ\P{p)) = ifpPii- 



(26) 



Substituting the meson state of Ea. (l23t into the above 
equation in the rest frame, and contracting the quark 
(antiquark) creation operators in the meson state with 
quark (antiquark ) annihilation operators in the quark 
field of the axial current, we can get the expression of 
the pseudoscalar decay constant 



fp = 



1 A— r dk\k\H,{k) 

2 27r-^ mp Jq 

(^, +m,)(^Q+mQ)- |fcp 
^EgEgiEg + mg){EQ + mg) ' 



(27) 



where Eq and Eg are the energy of the heavy and light 
quarks. To be consistent with the wave equation, here 
both heavy and light quarks are taken to be on-shell. We 
assume that the decays of the heavy meson can be ap- 
proximately described in terms of on-shell valence quarks, 
although the sum of the four-momenta of the valence 
quarks are not equal to that of the meson because of the 
existence of the color field within the hadron which can 
carry both energy and momentum. 

We would like to mention that the leptonic decay of 
pseudoscalar mesons was considered several decades ago 
in a different method by assuming the coupling of meson 
with quark- antiquark pair 25]. 



IV The QCD-inspired potential 

The potential of strong interaction between the heavy 
quark and light antiquark is taken as a combination of 
a Coulomb term and a linear confining term inspired by 
QCD 



11, 1 



V{r) 



4a^(r) 



b r 



(28) 



The first term is the Coulomb term, which is consistent 
with one-gluon-exchange contribution for short distance 
calculated in perturbative QCD. The second term is the 
linear-confinement term, which generates confinement in 
long distance. The third term is a phenomenological con- 
stant, which is needed to reproduce the correct masses for 
heavy-light meson system. 



The running coupling constant as{Q'^) in momentum 
space with Nf quark flavors at large values of Q^, calcu- 
lated in lowest-order QCD, is 



127r 



(33-27V/)ln(Q2/A^)' 



(29) 



This behavior of the strong coupling can be parameter- 
ized in a simpler form which can be conveniently trans- 
formed into the r-space 11] 



QV47? 



(30) 



where ai are free parameters chosen to fit the behav- 
ior of as{Q^) given by perturbative QCD (Eq. ip^ '). As 
q2 _^ ^jjg coupling as{Q^) diverges, which is believed 
to be a signal of confinement. However, as — A^, per- 
turbative QCD can not apply, the behavior of as{Q^) at 
small given in Eq. cannot be the exact prediction 
of QCD. One can make other choice for the behavior of 
the strong coupling at small momentum transfer. As in 
Ref. [llj, we assume that the coupling as saturate as a 
critical value q;"'*"^'*\ where q;^'''*"^'*' = J^i'^i- In. prac- 
tice, only several ai are needed to be non-zero, which can 
fit the behavior of as{Q^) well at perturbative region, de- 
viation only occurs at small Q^. 

The transformation of as(Q^) by using Eq. (pO)) in- 
stead of Eq. (PI is [Tit 



ai 



dx. 



(31) 



Fig IT] is the behavior of as{r) with the parameters ai — 
0.15, a2 = 0.15, as = 0.20, and 71 = 1/2, 72 = \/l0/2, 
72 = \/1000/2, which is relevant to the critical value 



, critical 



0.5. 



. 8 
£ 0.7 
a" 0.6 
. 5 
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. 1 



0.5 



r (fm) 



FIG. 1: The behavior of Qs(r), with the parameters cti = 
0.15, Q2 = 0.15, 03 = 0.20, and 71 = 1/2, 72 = v^/2, 
72 = \/T000/2. 



V Numerical result and discussion 

The parameters are selected by comparing the pre- 
dicted heavy meson mass with the experimental data. 
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The recently measured values of the decay constants of 
D and Dg mesons can give a further constraint on the 
parameters. The parameters which we finally obtain are 



5== 0.10GeV^ 
mt = 4.98 GeV, 
m, = 0.30 GeV, 

^critical ^ 



3= -0.19GeV^ 
TTic = 1.54 GeV, 
TO„ = mrf = 0.08 GeV, 



(32) 



and L = 10 fm, N = 100. 

The masses and decay constants of the B, Bg, D and 
Ds mesons calculated with the above parameters will 
be given in the following. The masses are given in Ta- 
ble m Here we do not consider the contribution of spin- 
dependent interactions in our calculation, it may give 
errors about 100 ^ 200MeV for the masses. Varying the 
parameters may also give errors to the numerical values. 
We estimate the combination of both the errors can be 
about 7% for B and B, mesons, and 10% for D and Ds 



TABLE I: Masses of pseudoscalar heavy mesons calculated 
by solving the wave equation, and the comparison with ex- 
perimental data. The data is quoted from the Particle Data 
Group :26|. 





niB 




mo 




this work 
(GeV) 


5.25 ±0.37 


5.34 ±0.37 


1.86 ±0.19 


1.96 ±0.20 


Exp. 

(MeV) 


5279.17 
±0.29 


5366.3 
±0.6 


1869.6 
±0.16 


1968.47 
±0.33 



The decay constants obtained are 

/s = 198 ± 14 MeV, /s^ = 237 ± 17 MeV, 

fn = 208 ± 21 MeV, fo, ^ 256 ± 26 MeV. (33) 

The comparison of the decay constants obtained in this 
work with experimental data are given in Table |lTl Both 
the masses and decay constants obtained in this work 
can be well consistent with experiment. For fs and /b^, 
there are still no precise measured values in experiments 
yet. Our prediction can be tested in experiment in the 
future. 

The Icptonic decay rates of B meson relevant to the 
decay constant /s obtained in this work are 

Br{B+ e+v^) = (1.11 ± 0.26) x 10"", (34) 
Br{B+ /i+z/^) = (4.7 ± 1.1) x 10"^ (35) 
Br{B+ ^T+Vr) = {l.l ±0.2) xlQ-'^, (36) 

where the errors are mainly caused by the uncertainties 
of the decay constant fs and the CKM matrix element 
Vub- The value of \Vub\ is quoted from PDG j2^ 



At present the branching ratio of Br{B~^ — )■ t^i^t) has 
been measured in experiment. The results still suffer 
from large uncertainties. The measured value of Belle 



collaboration is Br{B+ r+Vr) = (1.79+;;;^^+^|'J) x 



-1-0.56-1-0.46 N 
-0.49-0.5lJ 

lO""* [28|, while the values of BABAR collaboration are 
Br{B+ T+Ur) = (0.9 ± 0.6 ± 0.1) x 10"^ ^ 

-1-0.9 ' ' 



and 

combined re- 



X 

(1.8+^:^ ± 0.4 ± 0.2) X 10-4 [30|. The 
suit of BABAR collaboration is Br{B+ r+i^r) = 
(1.2 ± 0.4 ± 0.3 ± 0.2) X 10--* [31]. Considering the large 
uncertainties of the experimental results, our predictied 
branching ratio of the decay mode B~^ — >■ T~^i/r is cosis- 
tent with the experimental data. 

A super B factory will come into operation with the 
designed peak luminosity in excess of 10"^^ cm-^S"^ at 
the T(4s) resonance in the next half decades 3^. The 
integrated luminosity of 75 ab"^ would be collected in 
five years of data taking. Then the branching ratios of 
— > T^Vr and /i^f^ can be measured at the Super5 
factory with precisions of up to 4% and 5% respectively 
[s^ l. Taking the value of the CKM matrix element \Vub\ 
as input, the decay constant of /s can be obtained at 
super5. 

The decay constants are also compared with previous 
theoretical results in Table IIIII Our predictions for the 
decay constants are quite different from previous results 
calculated in the relativistic potential model. For fu and 
fo^ , our results are larger than that in Ref [14| , while our 
results for fs and fs^ are smaller than theirs. 45 1 



TABLE II: Decay constants of pseudoscalar heavy mesons 
calculated in this work, and the comparison with experimental 
data. All values are in units of MeV. 





fB 




Id 




this work 


198 ± 14 


237 ± 17 


208 ± 21 


256 ± 26 


Exp. [19, 27] 






205.8 ±8.5 ±2.5 


254.6 ±5.9 



The reduced wave functions in coordinate and momen- 
tum spaces (Eqs. ([8|) and (f20|) are depicted in Figs. [2] 
andO The wave function squared |u;(r)p is the possi- 
bility density distributed along the quark-antiquark dis- 
tance r. The curves in Fig. [2] show that the most proba- 
ble distribution of the quarks occurs at the distance 0.4 
fm between the quark and antiquark in both B and D 
mesons. The possibility density vanishes as the distance 
larger than 2 fm. The mean square root of the distance 
is about 0.5 ^ 0.7 fm. 

The numerical solution of the wave function in momen- 
tum space is given in Figl3l which shows that the peak 
of the momentum-distribution of the quarks in the heavy 
meson is at about 0.4 GeV. The reduced wave function 
can be fitted with the analytical form as suggested in 
Ref.fii 



\Vub\ = (3.93 ±0.36) X 10" 



ak 



(37) 
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TABLE III: The comparison of the decay constants calculated 
in this work with other theoretical results. All values are in 
units of MeV. 





/s 




fo 


Id, 


this work 


198 ± 14 


237 ± 17 


208 ± 21 


256 ± 26 


Ref.fl4]" 


230 ± 35 


245 ± 37 


182 ± 27 


199 ± 30 


Rcfjaa]" 


203 ± 23 


236 ± 30 


205 ± 20 


235 ± 24 


Ref.f34]'' 


206 ± 20 




195 ± 20 




Ref.f35]" 




— 


200 
230 


221 
270 


Ref.[36]" 






177 ± 21 


205 ± 22 


Rei.\37Y 


193 


195 


238 


241 


Ref.fSS]" 


196 ± 29 


216 ± 32 


230 ± 25 


248 ± 27 


Ref.[39]^ 


189 


218 


234 


268 


Ref.[40]'' 


210 ± 11.4 
±5.7 








Ref.[41j'' 






201 ±3 
±17 


249 ±3 
±16 


Ref.[42]'' 


195 ± 11 


243 ± 11 


207 ± 11 


249 ± 11 


Ref.[43]" 






207 ±4 


241 ± 3 


Ref.[44]'' 








248.0 ± 2.5 



a Relativistic potential model. 
b QCD sum rule, 
c Light front quark model. 
d Finite energy sum rules. 

e Quark model based on Bethe-Salpeter equation. 
/ Relativistic constituent Quark Model. 
g Derived from result of Lattice QCD. 
h Lattice QCD 




FIG. 2: The reduced wave function uo{r) in coordinate space, 
(a) is for B and Bs mesons. The solid curve is for the wave 
function of B meson, the dashed one is for B^. (b) is for D 
and Ds mesons. The solid curve is for D, and the dashed one 
for Ds. 



wave functions obtained here can be useful for studying 
heavy meson decays. 



where mn is the heavy meson mass. The factor 
AiTy/ mnOL^ is the normahzation factor due to the normal- 
ization condition in Eg. (1^51) . Note that the wave func- 
tion for B and/or D meson is 5'o(A) = (fo{k)/k. Our 
numerical solution gives a — 3.0 GeY-\ 2.6 GeV-\ 
3.4 GeV"\ and 3.2 GeV"^ for B, B,, D and mesons, 
respectively. 

With the constraint of the measured values of fo and 
fo^ considered, the wave functions obtained here can 
be more reliable than before, which should be useful in 
studying the decays of the B and D mesons. The appli- 
cation of the wave functions in studying the heavy meson 
decays deserves a separate work. 



VI Summary 

The wave functions and decay constants of B and D 
mesons are revisited in the relativistic potential model. 
The parameters in the potential model are further con- 
strained with the experimental values of fo and fo^ 
available. The masses and decay constants of the heavy 
mesons are obtained, which can be well consistent with 
the current experimental data. The wave functions both 
in coordinate and momentum spaces are obtained. The 
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